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At the 2004 AIAA/CEAS Aeroacoustic Conference, a breakthrough in acoustic 
microphone array technology was reported by the authors. A Deconvolution Approach for 
the Mapping of Acoustic Sources (DAMAS) was developed which decouples the array design 
and processing influence from the noise being measured, using a simple and robust 
algorithm. For several prior airframe noise studies, it was shown to permit an unambiguous 
and accurate determination of acoustic source position and strength. As a follow-on effort, 
this paper examines the technique for three-dimensional (3D) applications. First, the 
beamforming ability for arrays, of different size and design, to focus longitudinally and 
laterally is examined for a range of source positions and frequency. Advantage is found for 
larger array designs with higher density microphone distributions towards the center. After 
defining a 3D grid generalized with respect to the array’s beamforming characteristics, 
DAMAS is employed in simulated and experimental noise test cases. It is found that spatial 
resolution is much less sharp in the longitudinal direction in front of the array compared to 
side-to-side lateral resolution. 3D DAMAS becomes useful for sufficiently large arrays at 
sufficiently high frequency. But, such can be a challenge to computational capabilities, with 
regard to the required expanse and number of grid points. Also, larger arrays can strain 
basic physical modeling assumptions that DAMAS and all traditional array methodologies 
use. An important experimental result is that turbulent shear layers can negatively impact 
attainable beamforming resolution. Still, the usefulness of 3D DAMAS is demonstrated by 
the measurement of landing gear noise source distributions in a difficult hard-wall wind 
tunnel environment. 


Nomenclature 

a m shear layer refraction amplitude correction for e m 
A DAMAS matrix with A nn / components 

A nn / reciprocal influence of beamforming characteristics between grid points 

B f longitudinal beamwidth, R 2 — R { 

B' forward-defined longitudinal beamwidth, 2 (R 0 —R { ) 

B Wo lateral beamwidth of 3 dB down from beam peak maximum, = B of Ref. 1 
c 0 medium speed of sound 

CSM cross spectral matrix 

D diameter of array 

DR diagonal removal of G in array processing 
e steering vector for array to focus location 

e m component of e for microphone m 

f frequency 
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A f frequency bandwidth resolution of spectra 

FFT Fast Fourier Transform 

(j) array elevation angle 

G mm ' cross-spectrum between p m and p m/ 

G matrix (CSM) of cross-spectrum elements G mm ' 
y lateral resolution 

H height of chosen scanning plane 

i iteration number 

K acoustic wave number, 2 7tf t 

LADA Large Aperture Directional Array 

X longitudinal resolution 

m microphone identity number in array 

m 0 total number of microphones in array 

n grid point number on scanning plane(s) 

M wind tunnel test Mach number 

MEMS Micro-Electro-Mechanical System (microphone sensors) 

N total number of grid points over scanning plane(s) 
p m pressure time records from microphone m 
P m Fourier Transform of p m 

QFF Quiet Flow Facility 

Q n idealized P m for modeled source at n for quiescent acoustic medium 

r radial position on array face with respect to center 

r distance r m for m being the center c microphone 

X m retarded coordinate distance to m , T m c 0 

r 0 radius of array 

R distance of array from scanning plane 

R 0 distance of array to source location 

i?j distance of array to 3dB-down array response location forward of source 

R 2 distance of array to 3dB-down array response location aft of source 

SADA Small Aperture Directional Array 
STD standard or classical array processing 
T complex transpose (superscript), also block of time data 
TE trailing edge 

T m propagation time from grid point to microphone m 

W width of scanning plane 

Ax widthwise spacing of grid points 

x,y,z array coordinates 

x T ,y T ,z T tunnel coordinates 

X matrix of X n terms 

v, “noise source” at grid point n with levels defined at array, Q* Q n 

Ay heightwise spacing of grid points 

Y (< e ) output power response of the array at focus location 

Y matrix of Y n terms 

Y n Y ( e ) , when focused at grid point n 

A z longitudinal spacing of grid points 
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I. Introduction 

Traditional microphone array results represent noise sources that are convolved with array beamform response 
functions, which depend on array geometry, size, source positions, and frequency. The Deconvolution Approach for 
the Mapping of Acoustic Sources (DAMAS) method 1 represents a significant technology advance by removing 
these beamforming characteristics from output presentations. This can render an accurate and explicit definition of 
source strength and location. The authors in Ref. 1 dealt with DAMAS processing only for single two-dimensional 
(2D) scanning planes. The current study takes the first step at three-dimensional (3D) DAMAS processing. 

One can readily perform 3D processing with DAMAS simply by defining evaluation grid points on other than 
on a single plane. A primary issue that arises is one of attainable DAMAS resolution, which in turn depends on 
conventional beamforming resolution. For 2D lateral beamforming where a planar array focuses forward over a 
plane parallel to the face of the array, it was found 1 that DAMAS can spatially resolve as small as 5 to 20% of a 
lateral beamwidth (defined by 3 dB down array response). The smaller arrays have little or no focusing ability 
depth-wise (or longitudinally), so the “longitudinal beamwidth” may not be finite. The larger the array, the higher 
the ability is to focus depth-wise. However, too large an array would present a measurement dilemma in 
aeroacoustics, in that the array would enclose differing directivity regions of a source being measured, thus 
nullifying the measured results (since uniform radiation of any source being measured is assumed over the array 
traditional processing). In addition to size, array sensor distribution design also affects both lateral as well as 
longitudinal focusing characteristics. 

The present paper makes a first examination of 3D beamforming and corresponding DAMAS source 
distribution solutions, with regard to array design and size, as well as grid distribution. Classical beamforming is 
examined with respect to its ability to focus longitudinally and laterally as a function of array design and size with 
respect to source position and frequency. Three existing experimental designs are examined in some detail. Point 
source simulations and experimental data using a generalized grid are used to demonstrate attainable 3D resolution 
for beamforming and DAMAS. Array data obtained during aeroacoustic testing of a flap is used to examine 
localized and distributed sources. The role of test conditions such as the presence of turbulent shear layer on 
attainable spatial resolution is observed. Also, examined is a landing gear configuration in a hard wall wind tunnel 
environment with an array mounted on the sidewall. 

II. Basic Beamforming and DAMAS 

Figure 1 illustrates a generalized arrangement for a 
phased array distribution of m 0 microphones located outside 
a flow field containing an aeroacoustic model. A scanning 
plane of grid points is defined over the noise source region. 

Additional planes are shown that may be used when 3D 
distributions of sources are to be accounted for. The 
beamforming approach involves the generation of steering 
vectors associated with each microphone with respect to the 
chosen steering location. 

In Fig. 1, the steering location is designated as grid 
point n . The steering vector is 

e = COl [e, e 2 ••• e m J (1) 

where the component for each microphone m is 

e m =a m yexp{j27tfT m } ( 2 ) 

c 

The vector components serve to phase shift each 
microphone signal to allow constructive summing of 
contributions from the chosen locations. The steering vector contain terms that account for the mean amplitude a m 
and phase changes due to convected and refracted sound transmission through the shear layer to each microphone. 


additional scanning 
planes when 
3D imaging 



Phased Array 
of Microphones 


Scanning Plane 
n = 1 over Noise 

Source Region 


Figure 1. Illustration of open test configuration 
where the microphone array is placed outside the 
flow region containing the noise source scanning 
plane region. 
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T m is the time required to propagate from grid point n to microphone m . In retarded coordinates, \ m is the 
distance from microphone to the scanning location and r is that for the reference center microphone. 

The cross-spectrum matrix (CSM) is composed of cross spectra elements and is 1 


G = 



(3) 


where each term is the cross-spectrum between each combination of the array microphone using Fast Fourier 
Transformed (FFT) processed data. For classical or standard (STD) array beamforming, the output power spectrum 
(or response) of the array is obtained from 


Y(e)= 


\ T Ge 


m n 


(4) 


Here Y(e ) is a mean-pressure-squared per frequency 
bandwidth quantity. The division by the total number of 
array microphones M 0 -squared serves to reference levels 

to that of an equivalent single microphone measurement. 
Figure 2 illustrates a typical array result presentation over 
a plane through a source region showing dB level 
contours of Y(e ). Shading algorithms 1,2 3 in connection 

with Eq. (4) can be used over distributions of array 
microphones to modify and control the array beampattern 
beamwidth B defined as the width across the main 
response lobe over which the sensing level is within 3 dB 
from the peak level. Also, to improve dynamic range and 
remove extraneous noise, the diagonal of the cross- 
spectral matrix can be removed to obtain a special 
“Diagonal Removal (DR)” form 1,4 of Eq. (4) . 



Array 


D 

_L 


Figure 2. Geometric parameters of the array and source 
region scanning plane for DAMAS applications. 


DAMAS inverse problem definition and solution. As detailed in Ref. 1, the purpose is to pose the array 
problem such that one can back-out a unique distribution of sources that would fully explain the array output. 
Defining a pressure transform P m (for microphone m ) to be related to a modeled source located at position n in the 
source field, one gets 


P,, =Q n e (5) 

where for the modeled source, X n oc Q n Q n is the pressure-squared amplitude of the modeled noise at position n . 
From this, one obtains the modeled microphone array cross-spectral matrix CSM, G n , for a single source located at 
n . Upon summing sources from all n , the total modeled CSM is G mod . Then the modeled array output is 








( 6 ) 


By equating 7 (e) with processed Y(e ) from measured data, we have, when all sources over all n are 

considered statistical independent, 
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AX=Y 


with A being defined to have the components 


A ,= — 


[I 


(7) 

(8a) 


m n 


where the bracketed term is given by 


(<' ) < ( e i') e ‘ 


«)< 


«)■<_ 


(8b) 


Equation (7) represents a system of linear equations relating a spatial field of point locations, with beamformed 
array-output responses Y n , to equivalent source distributions X n at the same point locations. Reference 1 gives the 
details of the derivation above as well as the iterative solution method. 

2D Application Simulation. Figure 2 shows parameters found to be important in defining requirements for the 
DAMAS calculations. It was found 1 that a grid step size Ax (or Ay) satisfying 


0.05 < Ax/5(or Ay/B) < 0.2 


(9) 


is representative of the resolution range attainable 
the resolution as y ~ 0.155. Much finer resolution 
resolution would give less detail than needed and 
processing) in the DAMAS solutions. The spatial 
extent ratio is required to be at least W IB (and 
H / B) > 1, because it must be large enough that 
discrimination is apparent over the grid range in 
order that the grid points represent independent 
contributions to the DAMAS solution. With the 
variation over a distance B of only 3 dB, larger 
expanses are better. 

Figure 3 shows the example 1 of 2D 
beamforming and DAMAS for an ideal 
(synthetically generated) point source at R = 60” 
on a W = H = 50” plane. The array, SADA, is 
7.8” in diameter (its design will be detailed 
subsequently). With Ax = Ay=l” and beamwidth 
B ~ 1 foot, Ax/ B = .083. The top left frame of 
Fig. 3 shows a contour of array response in dB, 
when the array is focused at different grid 
locations over the plane, using standard (STD) 
processing of Eq. (4). DAMAS results are shown 
for one iteration (i= 1) of the DAMAS solver 
over the plane, for /=1000 iterations, and for 
z=5000 iterations. It is seen that progressively 
better, sharper DAMAS solutions for the point 
source is attained with more iterations. Sharper 
results are attainable for somewhat coarser grid 
spacing Ax. 


for a physical distribution of sources. Nominally one can denote 
would give more detail than physically feasible and much coarser 
perhaps cause aliasing error (in analogy with that of FFT signal 


dB 60 70 80 90 100 



_DAMAS, 1000 Iterations 


DAMAS, 5000 Iterations 


20 


S panwise location (in) 


Figure 3. Synthetic point source - SADA STD processing for 
50”x50” scanning plane that is 5’ away. Resolution Ax=1”. 
Frequency /= 10 kHz and Ax 75=0.083. 
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III. 3D Beamforming 


SAD A and LADA. The Small Aperture Directional Array (SAD A 2 ) of diameter D = 7.8” and the Large 
Aperture Directional Array (LADA 2 ) of diameter D = 34” are used in this section to illustrate 3D beamforming. 
(Note that the microphone patterns of these two arrays will be given in the following section.) Whereas Fig. 3 shows 
a lateral cut through a source at a distance R 0 = 60” from the center of and perpendicular to the face of the SAD A, 
Fig. 4 shows a longitudinal cut passing through both the source and the SAD A center. The levels are referenced to 
zero dB response when the array is focused to the point source. The 3 dB-down locations in front of the source 
toward the array, designated as Rj, and 3 dB-down to the sides of the source, whose distance between is B Wo , are 
reference measures of the beam pattern. In Fig 2 and previous discussions, the lateral B Wo is referred to as simply B. 
With the SAD A, the large ratio of Rq/D (=1.1) renders broad (large B Wo ) and elongated beam patterns, which do not 
“close” (i.e., do not form a closed beampattern) except at very large z distances (not shown in this scale) for the 
highest frequencies. In Fig. 5, corresponding results for the LADA (where R^/D =1.11) are shown. For all 
frequencies, the beam pattern is more tightly “focused”, with R 1 being closer to R 0 and the pattern less elongated. At 
6 kHz, the pattern is seen tighten, compared to that at 1 kHz. It is seen to “close”, or drop in level, at z values beyond 
the source location. There the 3 dB-down location is designated by R 2 . At higher frequencies of 10 and 40 kHz, the 
beam patterns are seen to tighten considerably (note the change of scales). Tables 1 and 2, to be subsequently 
discussed, include values for these reference length measures. 



Figure 4. SADA beamforming along a longitudinal plane 
with point source at R 0 = 60” (for f =1, 6, 10 and 40 kHz). 



Figure 5. LADA beamforming along a longitudinal plane 
with point source at R 0 =60” (for f =1, 6, 10 and 40 kHz). 


Beamwidth Parameters and Array Configurations. A sketch of a beam pattern defined by the 3 dB-down 
response curve is given in Fig. 6 for an array of diameter D, showing the lateral beamwidth B Wo at R 0 and locations 
R 7 and R 2 . A longitudinal beamwidth is defined as 


B,-R 2 -Ri 


and, also, an alternate definition is the “forward-defined” longitudinal beamwidth B' n 


( 10 ) 


6 

American Institute of Aeronautics and Astronautics 


B\ = 2(R 0 -R l ) 


( 11 ) 


The following analysis determines relations for these beamwidths as function of array type and source distance. This 
analysis is facilitated by viewing the array as a circular disk with a continuous sensing surface. This also permits an 
examination of array weighting (or microphone density) sensitivity with regard to an array’s ability to focus in 3D. 

For the circular array “disk” in Fig. 6, the Fourier transform 
F=Ffl ) of the resultant acoustic load response is related to the integral 
of the Fourier transform P(f,r, 0) of the acoustic pressure at radius r and 
polar angle 0over the disk, 

F= jj W(i\ d)P{f , r, 0)E(r, d)rdrdd ( 1 2) 

where W(r,0) is an array “sensitivity” weighting and E(f,r,0) is an 
“array steering” term. W(r,0) can be used to define a microphone 
density over the disk. E(f,r, 6) can be used to focus array response to 
specific locations in the sound field. The array disk response in terms 
of dB is 


dB = 10 Log(F*F) (13) 

The term F*F replaces Y ( e ) of Eq. (4) for the present continuous 

disk “array”. Equations (12) and (13) were evaluated for a number of 
array design weightings W(r) (microphone density variations) to 
determine the likelihood of the existence of optimum configurations 
for longitudinal focusing. The evaluations were carried out numerically using a mathematical software package and, 
alternately, using dense sensor distributions with standard array programming. 

Figure 7 shows sketches of these array disk weightings 
considered along with the LADA, SADA, and MEMS array (to be 
discussed). The Uniform distribution, for equal sensitivity across 
the array, is contrasted with the Conical distribution (in 
proportion to r/r 0 ) with zero sensitivity at the disk center r = 0 
and the greatest sensitivity at the edge at r = r 0 , the Inverse 
Conical (1 -r/r 0 ), as well as the inverse weightings (r/r 0 )" 7 and 
(r/r 0 )' 3/2 with sensitivities greatest approaching the center r = 0. In 
Fig. 8, array responses are shown, as represented by the 3 dB- 
down contours from the peak level at R 0 for the array designs. All 
diameters D were set to that of the LADA for comparison 
purposes, R 0 = 5’, and /= 10 kHz. In Fig. 8 (a) which shows the 
curves to the correct relative scale, it is seen that for all designs, 
the longitudinal beamwidths B } are large compared to their lateral 
beamwidths B Wo . This is disappointing to the extent it would be 
desirable to attain good and near-equal resolution of noise sources 
in all directions. Still there are important differences between the 
resolutions of the different designs, which may make 3D focusing 
practical or not depending on application. Differences are more 
easily seen in Fig. 8 (b) where the lateral coordinate is stretched. 
The array legends are listed in order of the ratio B l / B Wo that should be an indicator of how well the 3D focusing is 
for this f D , and R 0 case. The finite sensor number array designs do well in this comparison. The SADA has a 
relatively good B l / B Wo of 7.8. This is followed by the (r/r 0 )" 7 , (r/r 0 )' 3/2 , LADA, MEMS, Uniform, Conical, and 
Inverse Conical designs. 



INVERSE CONICAL INVERSE f 


90 ° 90 ° 90 ° 



LADA SADA MEMS 


Figure 7. Array weighting distribution 
functions considered over face of circular 
arrays. 



Figure 6. Sketch of a 3 dB-down beam 
pattern over a longitudinal plane for a 
point source at R 0 
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a. to SADA 


b. expanded scale 


Figure 8. Beam patterns of array response for arrays with differing weighting distributions (all are evaluated 
for R 0 ID = 1.765). 


Table 1 quantifies the lateral and longitudinal beamwidth values for the same array designs of Figs. 7 and 8, 
except SADA and the MEMS array are listed using their own diameters D. The values for SADA outer cluster, 
appropriate for the blended results presented in Ref. 1 at 10 kHz, are additionally presented. In Table 1, the lateral 
beamwidth B Wo is defined from the numerical calculations and the constant C 0 is subsequently defined using 


B 


w 0 


CqRq 

KD 


(14) 


Table 1. 3D Beamforming parameters for specific array designs for source at R 0 = 5’ and f = 10 kHz 


Array Disk Weighting 

D 

(ft.) 

Ri 

(ft.) 

r 2 

(ft.) 

Bi 

(ft.) 

Bi ' 

(ft.) 

-Cj 

B wo 

(ft.) 

Co 

B,/B wo 

Bi’IB wo 

UNIFORM 

2.83 

4.23 

6.05 

1.82 

1.54 

1.87 

0.15 

4.71 

12.3 

10.4 

CONICAL 

2.83 

4.18 

6.20 

2.02 

1.64 

2.08 

0.14 

4.45 

14.4 

11.7 

INVERSE CONICAL 

2.83 

4.10 

6.40 

2.30 

1.80 

2.33 

0.13 

4.10 

17.7 

13.8 

INVERSE r/ro 

2.83 

4.26 

6.00 

1.74 

1.48 

1.80 

0.18 

5.70 

9.7 

8.2 

INVERSE (r/ro) 3 2 

2.83 

4.20 

6.18 

1.98 

1.60 

2.05 

0.22 

6.90 

9.0 

7.3 

LADA - Unweighted 

2.83 

4.09 

6.38 

2.29 

1.82 

2.33 

0.22 

6.90 

10.6 

8.4 

SADA - Unweighted 

0.65 

1.01 

OO 

OO 

7.98 

2.30 

1.23 

8.98 

OO 

6.5 

SADA - Outer Cluster 

0.65 

1.04 

oo 

OO 

7.92 

2.30 

0.86 

6.28 

OO 

9.2 

MEMS - Unweighted 

2.46 

3.78 

6.91 

3.13 

2.44 

2.57 

0.24 

6.73 

13.0 

10.2 


Parameter Evaluation. The influence of array size, frequency, and source distance with respect to longitudinal 
beamforming are now examined. Referring to Fig. 6, for a point source at a distance R 0 , the pressure perceived over 

the face of the disk is P = P 0 exp (i(p 0 ), where the phase at radius r is (p^ — K(^Jr 2 +r 2 — R 0 ) with respect to that 
at the disk center where r = 0. This analysis assumes R 0 to be sufficiently larger than r to allow magnitude P 0 to be 
simply taken as uniform over the array. Specializing to the case of focusing the array to a point at distance R along 
the line through the array center and the source (thereby making the problem independent of 0), one can set, in Eq. 

(12), E — exp (iqj ) , where (p = K(^R 2 +r 2 — R) . The phase response over the disk at radius r is then 
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(15) 


(po ~ q> = K(^R\ + ? - R 0 ) - K(4r 2 +?-R) 

At R = R 0 , the phase would be zero over the whole array disk for maximum response for focusing at the source. It is 
hypothesized that in order for an array, focusing at R, to have a response (through Eq. (12)) that is a “set portion” of 
the maximum response, the phase (p 0 ~(p at any specific r would have a “set value”, which is dependent on the 
weighting W(r). Defining this “set portion” of response to be 3 dB down, and taking the specific r to be r 0 = D/2, 
and R = R x one has 


C, =K(^Rl +r 0 2 -R 0 )-K(^ 2 +r 0 2 -R t ) 

Equation (16) can be rearranged to give 


(16) 




(17) 


where the bracketed term [] is 



/ \ 

2 

2C, 

-l 

[2*0 | 

1 2R 0 
+ 1 + 0 

KD 

D J 

D 


(18) 


A corresponding term for R = R 2 is C 2 where C 2 is defined the same as C\ of Eqs. (1 6)-(l 8), but with subscript 2 
replacing 1 in both R x and C\. Note that C 2 = -C\. This relationship and the fact that the constants are virtually 
independent of frequency, array diameter, and R 0 , at sufficiently large RJD, were verified in numerical calculations. 
Values of C\ are given in Table 1 for each array design. 

As in Fig. 8, Table 1 considers only /= 10 kHz, although results for other frequencies can be determined from 
the above combination of equations. For the LADA and SAD A, Figs. (9) and (10) show values of B h B t \ B Wo , Bi 
/B Wo , and B{ !B Wo as functions of frequency f 




Figure 9. LADA lateral and longitudinal beamwidth 
parameters with respect to frequency. For R 0 ID = 
1.765 of Fig. 4. 


Figure 10. SADA lateral and longitudinal beamwidth 
parameters with respect to frequency. For R 0 ID = 7.69 of 
Fig. 5. 
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SADA 


It is seen that the 10 kHz cases of Fig. (8) and Table 1 show 
only intermediate relationships for the parameters that stabilize at 
the highest frequencies. For the LADA in Fig. (9), Bi is seen to 
approach a value somewhat less than For the SADA, B\ is off- 
scale. Figure (11) shows both Bi !B Wo and Bi IB Wo as a function of 
R 0 /D , for a number of KD values for the LADA sensor design. A 
corresponding figure for the SADA sensor design would be 
similar. The figure shows that Bi and Bi approach a set ratio with 
one another for high frequencies independent of the source 
distance from the array. 

It is seen that array design does affect 3D focusing 
characteristics as illustrated by Fig. 8. However, there is no 
obvious array sensor design that would allow very high-resolution 
3D focusing (small values of Bj !B Wo ) at practical values of RJD. 
In Fig. 11, the values of RJD are indicated for the LADA and the 
SADA in the QFF. This range is believed to represent a practical 
R 0 /D range in aeroacoustic testing, although for some applications, 
LADA is regarded as too large 4 . 



Figure 11. Ratio of beamwidth parameters with 
respect to R 0 ID for array KD values. Array 
microphone configuration is that of the LADA. 


IV. 3D Beamforming and DAMAS applications 


Grid definition and Simulations: Simulations for the LADA and SADA are performed to determine the 
degree to which 3D spatial source distributions can be resolved using DAMAS. Results are shown in Figs. 12-15. 



Figure 12. LADA longitudinal plane of 3D 

Beamforming output for synthetic point source at R„ Fi 9 ure 13 LADA 30 DAMAS solutions for cases of Fig. 

= 60”, for 1, 6, 10, and 40 kHz. Grid adapted to 12. Reduced ranges of z grid are shown, 

beamform parameters. 


The grid patterns are adapted to beamform parameters 
that are calculated by the preceding equations. The philosophy of the adapted grid pattern is consistent with that 
used in 2D applications 1 to meet the hypothesized needs for efficient DAMAS solutions. The grid spacing used is 
Kx/B Wo - AylB Wo = .167. For reference, the approximate lateral resolution is taken as this spacing, 

y ~ Ax = .167 B Wo which is consistent with that range given by Eq. (9). Because the resolution attainable in the 

/ 

longitudinal direction is the subject of study here, the ratio A z/B l was kept to a small .05 value in the vicinity of 
the source in order to be able to observe the resolution character. The width of the lower rectangular grid stack is 
kept at C Wo i? Wo , where C Wo is the number of lateral beamwidths that are chosen to enclose the source region 
containing R 0 . For the figures, C Wo = 3.5 is used. In this stack of 19 grid planes, the first plane starts at z=R=R 0 - 
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1.2(7?o-^i), others have Az=0.1(^ 0 -^i) spacing, and the last ends at ^=^ 0 +0.6(^ 0 -^i). Above the rectangular stack, 
the six expanding-width grid planes have widths equal to C Wo i? w , where B w is defined by Eq. (14) except R replaces 
R 0 . The positions for these grids are: first, R = Ro+0.8(R 0 -R\); second, ^ 0 + (Ro~Ri)l third, R 0 +l .5(R 0 -R\); forth, 
R 0 +2(R 0 -Ri); fifth, R 0 +4(R 0 -Ri); and last, 1.2^ 2 . However, if ^ 0 + 2(^ 0 -^i)^ R 2 , the second, third, and forth positions 
become, ^ 0 + A(R 0 -Ri)+A(R 2 - Ro), .5[i? 2 + Ro+ .4(i? 0 -^i) + -4(i? 2 - R 0 )], and R 2 , respectively. Also, if ^ 0 + 4(7? 0 -^i)^ R 2 , 
the fifth becomes 1.1 R 2 . Finally, the last grid plane position becomes 10^ 0 , if ^>10^0- 

For the FAD A, Fig. 12 shows a horizontal (longitudinal) plane cut of 3D beamforming levels for frequencies at 
1, 6, 10, and 40 kHz. Figure 13 shows the same for the corresponding DAMAS levels except the z scale is limited to 
100 inches and below in order to show details in the source region. For each frequency, the test grid pattern is square 
in x and y coordinates, so a vertical plane cut would give the same appearance. The levels are not contour plotted, 
but are represented by color for each block (the pertinent grid point for the block is located at the bottom left 
comer). The total number of DAMAS iterations employed over the 11025 grid points was i = 5000. More source 
dispersion was noted when i = 1000 iterations were used. With the exception of the 1 kHz results, DAMAS 
successfully recovers the source region near ^ 0 =60”. The peak at z = R 0 (where level is about 96.4 dB) and the 

adjoining grid points account for virtually all the noise source strength. Upon defining the resolution X as containing 

/ 

all but about 1 dB of the 100 dB source strength, it is seen that XI B l -0.15 (3 grid points) for frequencies of 6 kHz 

and above. For 1 kHz, substantial error occurs in the location and energy spread. In Table 2, the resolution results 

/ 

are summarized for these and several other frequencies. For FADA with ^ 0 =60”, it is seen that X IB l remains at 
about 0.15. The ratio of longitudinal to lateral resolution XI y is in the range of 7 to 8. In the table, the asterisk (*) 
indicates the entry is only suggestive or that there is no correct answer. This is seen to be the case for lower 
frequencies where Bj / B \ » 1 , where the longitudinal beamwidth “ceases to close”. 

For the SAD A, the beamform levels and DAMAS results are presented in Figs. 14 and 15. The grid pattern is 
defined using the same grid generation logic. Figure 15 limits the presented z scale to 200 inches. Except for the 
highest frequency, the source is not captured. Even then, the energy is substantially spread. The results are 
summarized in Table 2, where the results are consistent with the FADA results with respect to the condition 
B t / B \ » 1 , where there is failure of DAMAS to capture the source. 



Figure 14. SADA longitudinal plane of 3D Figure 15. SADA 3D DAMAS solutions for cases of Fig. 14. 
Beamforming output. Otherwise key as per Fig. 12. Reduced ranges of z grid are shown. 
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Table 2. 3D Beamforming and DAMAS point source simulation resolution results for the LADA and the SADA. R 0 =60”. 


Array 

Freq. 

(kHz) 

Ri 

(ft.) 

r 2 

(ft.) 

B, 

(in.) 

B,' 

(in.) 


B w0 

(in.) 

Az 

(in.) 

X 

(in.) 

AJB r 

A xl B wo 

Ax 

(in.) 

Aly 

LADA 

1 

1.33 

oo 

oo 

88.0 

oo 

26.0 

4.40 

35.2* 

0.40* 

0.167 

4.3 

8.1* 


3 

2.83 

17.16 

172.0 

52.2 

3.3 

8.7 

2.61 

10.4 

0.20 

0.167 

1.5 

7.2 


6 

3.64 

7.80 

49.9 

32.7 

1.5 

4.3 

1.63 

4.9 

0.15 

0.167 

0.7 

6.8 


10 

4.09 

6.38 

27.5 

21.8 

1.3 

2.6 

1.09 

3.3 

0.15 

0.167 

0.4 

7.5 


20 

4.50 

5.61 

13.3 

11.9 

1.1 

1.3 

0.60 

1.8 

0.15 

0.167 

0.2 

8.2 


40 

4.74 

5.29 

6.6 

6.3 

1.1 

0.7 

0.31 

0.9 

0.15 

0.167 

0.1 

8.6 

SADA 

1 

* 

oo 

oo 

* 

oo 

147.7 

6.10 

* 

* 

0.167 

24.6 

* 


3 

0.29 

oo 

oo 

113.2 

oo 

49.2 

5.66 

* 

* 

0.167 

8.2 

* 


6 

0.63 

oo 

oo 

104.9 

oo 

24.6 

5.24 

* 

* 

0.167 

4.1 

* 


10 

1.00 

oo 

oo 

96.0 

oo 

14.8 

4.80 

33.6* 

0.55* 

0.167 

2.5 

21 .5’ 


20 

1.69 

oo 

oo 

79.6 

oo 

7.4 

3.98 

23.9* 

0.35* 

0.167 

1.2 

22.6 5 


40 

2.53 

175. 

2070. 

59.3 

35 

3.7 

2.97 

14.8* 

0.25* 

0.167 

0.6 

24. T 


Experimental Cases 

In this section, experimental data are used to validate the resolution results from the simulations and to apply 
3D beamforming and DAMAS for several wind tunnel aeroacoustic test setups. Data used are from LADA 
measurements for a calibration “point source” and from a flap/airfoil model test conducted in the Quiet Flow 
Facility (QFF) at NASA LaRC. Also used are data acquired with the NASA LaRC MEMS sensor array during 
landing gear testing in the Virginia Tech (VT) Stability Wind Tunnel. For these cases, unlike the pure tone 
simulation cases, the results are presented in terms of one-third octave levels. This requires consideration of 
bandwidth error issues. 

Spectral bandwidth. For the array testing 1 done in the QFF, the data has been processed by Fast Fourier 
Transform (FFT) methods rendering a frequency resolution of A f = 17.44 Hz = l/T, where T is the individual data 
ensemble block size, in units of seconds. For the VT test data, the A f used was 24.4 Hz. For each cross-spectral 
component of the CSM, Eq. (3), there is a center frequency and phase associated with each A f band. With A f 
being sufficiently small, there is no error issue with the fact that Af ^ 0 . But with A f being small, repeated 
applications of DAMAS at each Af are required, and then summed, to produce single one-third octave DAMAS 
source maps. This is computationally time consuming, so it is desired to increase the bandwidths Af so that less 
summations are required to form one-third octave results. The question arises of the error associated with the chosen 
sizes of Af and how the bands are formed. 

Considered here are cross-spectral errors associated with time delay differences between microphones. First, 
there is a negative error due to processing when there is a time delay T between microphones m k and m l . The 
delay T is dependent on source location, the relative microphone locations, the test median speed-of-sound, and any 
required shear layer corrections. The FFT processing produces an “estimated” cross spectrum 5 


(G v ,\ st =V-pG m 


(19) 


where is the “actual” cross spectrum. The error is independent of frequency. The error is highest between 

extreme microphones of the array when focusing to edges of extended scanning planes. The effect is to effectively 
attenuate the contributions of the outer edges of the array - in contrast to the full coherence that is assumed in matrix 
A of Eq. (7). Equation (19) can be used to guide one’s choice of reprocessing with different T block sizes in order 
to maximize A f, and still keep error within chosen bounds. To keep the error below 0.2 dB (a stringent choice) in 


12 

American Institute of Aeronautics and Astronautics 



any of the CSM components, one can reprocess to attain A f -220 Hz for the SADA setup geometry in the QFF with 
a scanning plan at R 0 = 60” and 50” wide. For the LADA, A f -50 Hz for the same scanning plane and Af ~ 100 Hz 
for a one half size scanning plane (more practical anyway for the LADA in terms of number of required grid points, 
considering Ax requirements with respect to LADA’s smaller B Wo values). 

An alternative to reprocessing the cross spectra from the original time signals in order to increase A/ is to 
combine existing Af bands vectorially. There is a bias error inherent in doing this, due to variation in phase because 
of the same time delay T considered above. An “estimated” cross spectrum with new Af created by summing 
bands, would be 


1 | 

) est A r J G mm ff 

A/ fx 


sin (ttA/t) 

G = 

sin (7TT/T) 

nkfT 

m k m i 

_ 7TT/T _ 


( 20 ) 


-<b=90° 


>= 106 ° 


LADA 


where the form for the cross spectrum is taken as G m ^ — | exp(z'2 7Jff) . Here, Af — f — f and T is the block 

size it would have taken to obtain Af in normal processing. The bracketed term of Eq. (20) is substantially closer to 
the value 1 than that of Eq. (19) at small T IT . For the same level of CSM error, the corresponding Af values 
possible for the SADA and the LADA with Eq. (20), are now 3.4 times larger now at A f ~ 750, 170, and 340 Hz 

compared to the A f~ 220, 50, and 100 Hz, respectively, 
given in the last paragraph for Eq. (19). 

In the following applications for the QFF and VT tests, 
Af were increased by summing bands in accordance with 
Eq. (20), although a more conservative approach was 
actually taken by using only about a half of the Af 
permitted as guided by the examples given above. In a 
separate matter, it should be noted that the permissible Af 
values, found using Eqs. (19) and (20), are independent of 
frequency f One could use other criterion for analysis 
bandwidths Af such as those for 1/1 2 th octave bands 
(Af-. 058//) as long as permissible chosen values of Af 
are not exceeded at each DAMAS calculation. 

QFF “Point” Source Test. A sketch of the LADA 
positioned outside of the flow field in the QFF is shown in 
Fig. 16 (a) for the model flap and main airfoil element test 
set-up. The scan planes shown are for the subsequent study 
of flap edge and flap cove noise. The airfoil is at 16° angle- 
of-attack to the vertical. A similar experimental set-up using 
the SADA was reported in Ref. 1 . The LADA is positioned 
at an angle of (f) = 106°, referenced to the vertical line 
through the main element trailing edge (TE) and at a 
distance of R 0 = 5 feet. The use of the LADA is detailed in 
Refs. 2. Its microphone arrangement is shown in Fig. 7. 
Figure 16 (b) shows the calibrator source arrangement 
located at the TE next to the flap edge. The source is an 
open end of a one-inch diameter tube emitting broadband 
noise. 

Figures 17 and 18 show, respectively, the 3D 
beamforming and DAMAS for the calibrator source output 
with no tunnel flow (i.e., M = 0) and for the same grid 
pattern as the simulation of Figs. 12 and 13. As mentioned, 
the presentations are for one-third octave frequency bands, 
obtained by combining CSM results from a number of 



a. Noise path from flap edge to the LADA. Scan 
planes are for flap noise study. 


Calibrator Source 



Model 


b. Open end of calibration source positioned next to 
flap edge. 


Figure 16. LADA test set-up in QFF for flap noise 
test and calibration. 
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individual 17. 44 Hz bands to obtain Af=l22 Hz, used for evaluations of DAMAS. The DAMAS results were then 
summed to obtain the one-third octave results. The experimental results displayed in Figs. 17 and 18 show there is 
an approximate 0.3” off-set in the assumed source center position and some broadening (about 0.4”) due to the 
finiteness (1” diameter) of the source. But there is also some additional longitudinal smearing of about 0.05 to 0.1 
(1 to 2 blocks) over that found for the synthetic point source results displayed in Figs. 12 and 13. The longitudinal 
resolution found here is about 6”, 5”, and 1” for 6, 10, and 40 kHz, respectively. 



Figure 17. LADA 3D beamforming output for Figure 18. LADA 3D DAMAS results corresponding 

experimental calibration source of Fig. 16 (b) to Fig. 17 

presented over same grid and format of Fig. 12. 


QFF Flap Test. The flap edge test configuration for the LADA in the QFF is shown in Fig. 16 (a). The test 
case considered is for a 29° flap angle setting and M= 0.11. For a similar test case (see Ref. 1) in which the SADA 
was used, it was shown there was strong localized flap edge noise and distributed flap cove noise. This case is used 
here for a 3D DAMAS evaluation for multiple sources of different distributions. Of course, it is already known that 
both sources can be regarded as being located along one scanning plane that is aligned with the airfoil main element 
chordline. This scanning plane is designated in Fig. 16 (a) by R 0 (= 60”). For this scanning plane, Fig. 19 shows 2D 
beamforming contours and DAMAS results over the flap region using Ax = Ay = 0.25” for the LADA with f 1/3 = 20 

kHz (A xl B Wo = .19). For this M ^ 0 tunnel case, Fig. 19 shows results using Diagonal Removal (DR) processing for 
the first time in this report. The region of the contour plane that is presented as blank actually contains negative 
pressure-squared beamforming output from the DR processing. DR DAMAS processing was shown in Ref. 1 to give 
less noisy results by reducing the contributions of extraneous microphone wind self noise from the calculations. 
(Still, this result shows substantial extraneous “noise” over the scanning plane. This is the result of the array peering 
through the shear layer turbulence 4 , not microphone self noise. This “noise” is more severe for larger arrays. Note 
that if only 2D results were needed for this measurement exercise then the SADA or a sub-array of the LADA would 
be preferred.) The strong localized source region at the flap edge and the distributed source along the flap and main 
airfoil element cove are captured by 2D DAMAS. The total flap edge noise and flap cove noise levels are 44.5 and 
42.1 dB, respectively, (These were found for the flap edge by summing over x T = 39” to 44” and y T = ± 1” and for 
the cove by summing over x T = 38” to 41” andy T = 1” to 18”.) Note that by using LADA sub-arrays of various sizes, 
DAMAS determines the same total flap edge and flap cove noise levels over the same areas within 1 dB. For the 
flap edge noise, the above agrees with the 44.6 dB found for the SADA in Ref. 1 for the similar test case mentioned. 
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In order to perform 3D calculations, a target limit of 
20,000 grid points was used, which required a reduction of 
the scanning plane size from that of Fig. 19 and a reduction 
of the number of planes compared to Figs. 17 and 18. Figure 
16 (a) shows the scanning planes at R= 54” (=Ri), 57”, 60” 
(=5 0 ), 63”, 66”, and the last plane at 81” (=1.2 R 2 ), which 
partially follow the grid pattern logic previously used. The 
planes are normal to a line connecting the array center to the 
trailing edge (TE) and flap edge junction. In the array 
coordinate y direction shown, the grid used C Wo = 3.5. The 
grid is shown for the array coordinate y (= 71 ) spanwise 
dimension. Figure 20 shows beamforming results over the 
six scanning planes. The result for the R 0 plane is the same 
as that of Fig. 19. Although the levels are the highest at 5 0 , 
the levels on planes on either side are very close to these. 
Any de- focusing away from R 0 is barely apparent visually, 
until the last plane at 1.2 R 2 . The 3D DAMAS results of Fig. 
21, show a somewhat stronger distinction between the 
scanning planes, although clearly the sources are not 
isolated to the 5 0 plane. (Ideally, the test model sources in 
the R 0 plane would match those of Fig. 19 and not be 
included in the other planes.) The noise levels, summed in a 
manner consistent with that done for Fig. 19 and for the 
planes with increasing R values are 30.4, 37.8, 42.0, 24.5, 
32.0, and 24.2 dB for the flap edge, and 34.5, 34.6, 37.6, 
34.2, 34.2, and 29.0 dB for the flap cove. Compare these to 
the 2D totals of 44.5 and 42.1 dB for the flap edge and flap 
cove, respectively. 



S panwise location (in) 


B 


dB 

40 

37 

34 

31 

28 

15 

<0 


dB 


fl 


35 

25 

15 

5 

-5 


Figure 19. 2D beamforming output and DAMAS 
for flap noise test. DR processing with 

f m = 20kHz for Av = 0.25 


The results are good in that the highest levels are located on the R 0 plane. However, the results are problematic 
and important in that they indicate that the resolution scale A is larger than anticipated from the preceding 
simulations and calibrations. This can directly impact the usefulness of 3D measurements, at least for this test set- 
up. With the definition that A is the longitudinal scale that accounts for all but about 1 dB of source level (assumed 
to be found by the 2D results), one can find A by combining the pressure-squared values of the R 0 and adjacent 
planes. For the localized flap edge noise source, one obtains A ~ 0.255/, but for the distributed flap cove noise 
source, A ~ 0.45/. Here 5/ is calculated by the forgoing equations. The loss of resolution from A ~ 0.155/ is viewed 
by the authors as due to an effective increase in beamwidth 5/ due to the presence of the open-jet turbulent 
shearlayer. As was observed for Fig. 20, the beamformed pattern for the five adjacent planes are very similar with 
one another with little de-focusing apparent as one moves away from R 0 . More difference should be seen across the 
one beamwidth distance, B h that is spanned by the five planes. It is not an issue with DAMAS per se - but with 
beamforming. As discussed to some length in Ref. 4, sound produced by the model is transmitted through the 
turbulent shear layer on its way to the microphones. The ray paths are scattered by turbulence with respect to the 
mean refracted ray paths, thereby randomizing the noise received by each microphone. The importance of this 
scatter is more severe for the larger array sizes, because different microphones receive rays of different scatter paths 
thereby reducing coherence over the face of the array. This serves to spatially blur the source. This has not been 
considered a general problem for lateral beamforming for large arrays where resolution is already high. And no 
energy is lost in spatial integration 4 if one can increase spatial expanse. However, for longitudinal beamforming, 
resolution blurring is an issue if intrinsic beamwidth is already considered too large to define the source field 
properly. 
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3D beamforming output for flap noise test. DR processing with f m = 20 kHz. Ax = 0.25 ”. Five stacked scan 
one distant background scan plane. 
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Figure 21. 3D DAMAS results corresponding to Fig. 20. 
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VT Landing Gear Test. An aeroacoustic noise test 6 was 
conducted for a 26%-scale, high-fidelity Boeing 777 main 
landing gear (LG) model in the Stability Wind Tunnel at the 
Virginia Polytechnic Institute and State University (Virginia 
Tech, VT). The model was previously tested at VT 7 and at 
NASA Ames Research Center 8 . A portion of the VT 
measurements was conducted by NASA LaRC with the 
primary purpose of demonstrating the use of an array using 
Micro-Electro-Mechanical System (MEMS) sensors 6 . 


The 128 MEMS microphone array is shown mounted in 
the sidewall of the 6 by 6 foot test section mounted beside the 
LG in Fig. 21. A ceiling mount was also employed in the test. 
The array pattern is shown in Fig. 7. The array is recessed 
from the tunnel wall by 1” and is covered by a tautly drawn 
Kevlar material to maintain wall integrity and to prevent 
boundary layer turbulence from producing pressure fluctuations 
directly over the microphones. This reduced noise floor levels 
substantially. Because of the edge blockage about the array cut- 
out, the outer microphone ring is not employed making the 
array effectively 32” in diameter with an active xx 
microphones. Figure 22 shows a top-view of the LG image and 
the side-wall mounted MEMS array with locations indicated for 
scanning plane cuts used in the 3D beamforming and DAM AS 
calculations. As was done for the preceding flap calculations, to 
keep the number of grid points down to approximately 20,000, 
a reduced number of grid planes in the source region is used. 
But in contrast no distant grid planes are used. With the array at 
z = -1”, the stacked planes are atz= 19”, 24”, 29”, 33.5”, 38”, 
and 43”. 


Figure 21. MEMS array in side mounted with 
26% scale 777 landing gear test set-up in the 
Virginia Tech Stability Wind Tunnel. 
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Figures 23 and 24 show beamforming and DAMAS results 
from the side-wall mounted array for the LG with the wheel 
truck angle at 13° and the tunnel Mach number at M = 0.17. 

The one-third octave presentations are for f m = 6.3 kHz. With 
the array being relatively large and close to the distributed LG 
source region, the parameters cited as important for 3D DAMAS calculations are defined in a mean sense. Strictly 
speaking, because sources exist over each plane, each can be regarded as defining a R 0 position. However, upon 
choosing as a reference source distance R 0 ~ 1”+ 33.5”, one has R 0 /D ~ 1.08. Using preceding equations, one has 
Rj~ 26.9”, R 2 ~ 46.8”, Bi~ 20”, Bi’~ 15.2”, so the resolution is calculated as 3”. Nominally B Wo ~ 2.3”, so with 
Ax= .85”, one has A x/ B Wo ~ .37 (which is outside of y/B Wo range between .05 and .2). For the W= 50” and H= 
50”, chosen to cover the spatial range of the LG. The planes are large with respect to the lateral beamwidth, W/B Wo 
= C Wo ~ 22. The space between planes in the stack is Az= 4.5” and 5”, which exceeds the reference resolution X 
given above. (No accounting is attempted for any effective increases in the value for B t and subsequently X that may 
relate to the turbulent boundary layer scattering effects that was found present for the flap noise test application.) 
The present grid is seen to be a significant compromise that exceeds or are on the border of acceptable values of 
parameters important to 3D DAMAS application. For this reason, the application is considered preliminary. 


Kevlar Sheet f Recessed Array 


Figure 22. Stacked scan planes in 3D 
Beamforming and DAMAS for landing gear of 
Fig. 21 


In Fig. 23, the beamform results for the different cut planes show a substantial concentrated noise region toward 
the strut braces. In the planes farther from the array, sources more toward the wheel truck become apparent. Note 
that the steering vectors, Eq. (2), for the microphones accounts for shear layer refraction at the wall/tunnel flow 
interface at z = 0 and the mean tunnel flow convection speed. The effects of the Kevlar sheet at the shear layer 
interface has not been taken into account in the steering vectors. 
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Figure 23. 3D Beamforming results for side-mount MEMS array over 6 planes about the LG. 

The corresponding DAMAS results in Fig. 24 separate and quantify sources between the chosen planes. The 
first plane from the array at z = 19” cuts through the outer portion of the strut braces. Noise contributions from the 
braces appear to stand out, although there are sources assigned to the upper and center part of the plane that have to 
be regarded as being extraneously positioned. The next plane is at z = 24”, just in front of the wheels that are nearest 
the array. The plane cuts about midway through the strut brace assemblies, which seem to show strong distributed 
noise sources. For the third plane at z = 29”, just behind the near wheels, there appear to be contributions from the 
upper part of the strut braces and cross members. For the center plane for the truck assembly at z = 33.5”, and over 
planes just in front of and behind the far wheels at z = 38” and 43”, there seems to be contributions from the front 
brake and cable harness areas about the wheel truck, as well as the door panel and its mounts on the other side of the 
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stmt. It is remembered that DAMAS assigns dB level to sources with respect to the level perceived by the array. 
Therefore, for a source that is farther from the array to give the same level, as one that is closer, it would have a 
higher “source strength”. 
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Figure 24. 3D DAMAS results for side-mount MEMS array corresponding to Fig. 23. 


With respect to Fig. 24, there are sources distributed over the borders of the scanning planes, particularly along 
the bottom below the floor reference lines on the first and last scanning planes. This effect is discussed in Ref. 1, 
where it is mentioned that DAMAS apparently assigns “energy” or sources to borders when beamforming results 
suggest all sources are not within the defined scanning plane(s). Note that the floor region is a strong source of 
reflection that is not included in the particular grid planes chosen. Not included, also, is distant scanning planes 
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which potentially could be used to account for reflected / reverberant noise due to the tunnel walls. Some 
preliminary calculations indicated it would be counterproductive or wasted for this array and source orientation, 
which showed no identifiable strong reflection directions except for the afore-mentioned floor. The present 
application of Fig. 24 is not considered optimal. 


V. Conclusions 

Previously, beamforming and DAM AS post processing was dealt with only for single 2D scanning planes. The 
present paper examines the techniques for 3D applications. The basic beamforming ability of arrays, with regard to 
size and design, to focus longitudinally and laterally is evaluated as a function of frequency and distance. For all 
practical designs and measurement situations considered, it was found that the longitudinal beamwidths Bi are large 
compared to their lateral beamwidths B Wo . At 10 kHz, for a particular 34” diameter array at 60” from a source, Bi 
-27.5” and B Wo ~ 2.6”, so BJ B Wo ~ 10.6. This restricts one’s ability to attain good resolution in all directions. Still, 
array design does have an impact on measurement practicality. It was found that some existing arrays rate relatively 
well in this regard. Formulations are developed to predict attainable beamforming spatial resolution as a function of 
design. These were validated and incorporated in an evaluation of the ability of DAMAS to extract the source 
location and strength from beamformed results. It was found that the attainable spatial resolution X (containing all 
source energy within 1 dB) was about 15% of the longitudinal beamwidth B h This percentage is about the same for 
lateral resolution, for 2D and 3D, with respect to lateral beamwidth B Wo . An experimental calibration, which did not 
involve tunnel flow, validated this but found some additional spatial smearing to attain X of about 20% of the 
theoretical B h 

Experimental data from two different wind tunnel tests were examined in this study. A flap noise test was used 
to evaluate 3D DAMAS for multiple sources of different type distributions. The model was in an open-jet tunnel 
flow and the array was placed outside. It was found that the attainable resolution X was less that than ideal at about 
25% of the calculated B\ for a strong localized flap edge source. For a more broadly distributed flap cove noise, X 
was about 40% of B h This is believed to be the result of an effective beamforming related increase in B\ due to the 
presence of turbulence, and subsequent noise ray path scattering, in the open-jet shear layer. This is viewed as 
problem that limits somewhat the potential routine utility of 3D beamforming and subsequent DAMAS application. 
The effect would also be found in closed wall wind tunnels where arrays are mounted a sidewall with an associated 
turbulent boundary layer. 

The other experimental application is that of a landing gear in a closed wall windtunnel. In this case the array is 
relatively large and close with respect to distributed landing gear noise sources. This is not considered an optimum 
application, as there were some compromises made with regard to grid coverage and DAMAS application 
parameters. Still, the DAMAS results did successfully demonstrate the ability to separate and quantify sources over 
regions of the landing gear model. 
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